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The  optical  and  Zeeman  spectra  were  reported  for  Yb3  +  in  nine 
single  crystal  scheelite  hosts  in  a  previous  report. 
Identification  of  the  Si»  spectral  lines  allowed  the  energy  level 
schemes  for  Yb3+  in  the  nine  scheelites  to  be  determined.  From 
this,  best-fit  sets  of  crystal-field  parameters  were  calculated 
to  fit  the  expansion  of  the  crystal-field  potential,  Vc 
■  EB?'d5.  The  dependence  of  the  BjJ  on  the  unit  cell  volume  was 
examinee!,  and  it  was  found  that  the  strength  of  the  crystal  field 


1473  EDITION  OP  I  NOV  EE  IE  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OP  THIS  PAOE  fBRpii  Dpfc  Bntorod) 


A:  c. 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THU  FAPEOH—  D*m  »ll«m) 


20.  ABSTRACT  (Cont'd) 

as  measured  by  the  crystal-field  parameters  and  the  resulting 
Stark  splitting  decreased  when  the  unit  cell  volume  increased. 
It  was  found  also  that  the  dependence  of  the  Stark  splitting  was 
more  sensitive  to  the  lower  order  parameters,  B® ,  B®  ,  and  B[j  , 
than  to  the  higher  order  ones,  B®  and  Bf1* ,  by  almost  an  order  of 
magnitude. 

Crystal-field  parameters  were  also  calculated  from  first 
principles  by  using  the  conventional  electrostatic  point  charge 
model,  and  the  results  of  these  calculations  were  compared  with 
the  empirically  determined  parameters.  As  had  been  found 
previously,  agreement  was  poor  with  the  discrepancies  between 
experimental  fit  and  calculation  an  order  of  magnitude  or 
greater.  Various  phenomenological  corrections  to  the  point 
charge  model  were  investigated  including  (1)  adjustment  of  the 
ligand  oxygen  effective  charge  to  account  for  the  covalent  nature 
of  the  heavy  metal  tetrahedron,  (2)  application  of  an  effective 
dielectric  constant  depending  exponentially  on  distance  from  the 
impurity  ion  to  take  into  account  the  effects  of  lattice 
polarization,  (3)  expansion  of  the  4f  radial  wave  function 
(nephelauxetic  effect),  and  (4)  reduction  of  the  nearest  neighbor 
(ligand)  distance  to  account  for  the  local  lattice  distortion  due 
to  the  size  discrepancy  between  the  Yb3  +  ion  and  the  normal 
cation.  Crystal-field  parameters  were  calculated  by  adjusting 
various  parameters  to  account  for  these  corrections.  Good 
agreement  with  experimental  results  was  obtained,  the  discrepancy 
between  experiment  and  theory  being  reduced  to  about  8  percent. 
The  adjustment  of  the  ligand  oxygen  charges  and  the  expansion  of 
the  radial  wave  functions  proved  to  be  the  most  important 
modifications  tested  in  obtaining  the  fit. 
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1.  INTRODUCTION 

\ 

V  *  / 

In  a  previous  report*  (hereafter  referred  to  as  I),  we  presented 

the  results  of  a  series  of  measurements  of  the  optical  and  Zeeman  spec¬ 
tra  of  trivalent  ytterbium  (Yb3+)  doped  into  the  following  members  of 
the  scheelite  family  of  crystalline  hosts:  cadmium  molybdate  (CdMooV), 

p  p 

calcium  tungstate  ( CaW04 ) ,  calcium  molybdate  (CaMo04),  strontium  tung- 

p  n  p 

state  ( SrW04 ) ,  strontium  molybdate  ( SrMo04) ,  lead  tungstate  (PbWC>4) , 

r  ?  p 

lead  sK>lybdate  (FbMo04),  barium  tungstate  ( BaW04 ) ,  and  lithium  yttrium 
fluoride  (LiYF4).  Measurements  were  made  at  temperatures  varying  from 
below  the  lambda  point  of  liquid  helium  up  to  room  temperature  on 
samples  with  impurity  ion  concentrations  varying  from  0.05  to  4.0  per¬ 
cent  Yb.  The  results  were  the  identification  of  the  electronic  transi¬ 
tions,  both  absorption  and  fluorescence,  and  the  measurement  of  the  g 
factors  of  the  lowest  J  =  5/2  state.  These  results  are  given  in  tables 
3  and  4  of  I. 

The  objective  of  this  report  is  to  take  these  experimental  results 
and  from  them  derive  a  self-consistent  set  of  crystal-field 
parameters.  Following  this  calculation,  using  the  conventional  approach 
of  the  electrostatic  point  charge  model,  we  calculate  a  similar  set  of 
crystal-field  parameters  to  compare  with  the  experimentally  determined 
ones . „  Initially,  the  agreement  is  poor,  as  might  be  expected  from  the 
results  of  previous  works2”5  on  the  subject.  We  then  discuss  modifica¬ 
tions  to  the  point  charge  model  that  take  into  account  phenomena  that 
may  be  occurring  in  the  lattice.  These  modifications  include  (1) 
adjustment  of  the  ligand  charges  to  account  for  the  covalent  nature  of 
the  heavy  metal  tetrahedra,  (2)  expansion  of  the  4f  radial  wave  functon 
(nephelauxetic  effect),  and  (3)  reduction  of  the  ligand  distance  to 


*See  numbered  references  in  Literature  Cited  Section,  pp.  34,35. 
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account  for  local  distortion  in  the  lattice  due  to  the  size  discrepancy 
between  the  Yb3+  ion  and  the  normal  cation.  These  modif icatons ,  partic¬ 
ularly  accounting  for  the  nephelauxetic  effect,  improve  the  fit  to 
experiment. 

2.  DETERMINATION  OF  EXPERIMENTAL  CRYSTAL-FIELD  PARAMETERS 

The  trivalent  ytterbium  ion  has  a  4f13  configuration,  which  is  split 
by  the  spin-orbit  interaction  into  two  states,  2Fg^2  and  2F7^2,  the 
latter  having  the  Hund's  rule  ground  state.  These  two  states  are  then 
split  into  seven  levels,  each  doubly  degenerate  (Kramer's  degeneracy), 
by  the  electric  field  of  the  various  ions  that  make  up  the  lattice  in 
which  the  Yb3+  is  imbedded.  Group  theory  predicts  that,  for  an  Si* 
crystalline  field,  the  composition  of  the  two  spin-orbit  states  will  be 
2r5^6  +  2 T7f8  for  the  J  *  7/2  level  and  2r5f6  +  r7j8  for  the  J  *  5/2 
level,  thus  giving  a  total  of  seven  states. 

The  Hamiltonian  representing  the  crystal  field  can  be  written  as  a 
multipole  expansion  of  an  electrostatic  field  expressed  in  terms  of 
spherical  tensors.  In  the  notation  of  Wybourne,6  this  becomes 

vo  *  I  ■£*<£♦  »($!)  .  (1) 

n,m 

where  the  coefficients  of  the  tensors  are  called  crystal-field  parame¬ 
ters.  The  second  term  represents  the  spin-orbit  interaction  and,  to  a 
first  approximation,  could  be  ignored  since  the  average  Stark  splitting 
is  300  to  500  cm"1  compared  with  ~10,000  cm"1  spin-orbit  splitting. 
Thus,  J-mixing  has  little  effect  on  the  energy  levels.  However,  we  do 
include  the  spin-orbit  term  since  it  was  found  that  J-mixing  has  a 
profound  effect  on  the  g  factors. 


The  sum  in  equation  (1)  can  be  written  explicitly  as  follows . 
First,  for  a  lanthanide  rare  earth  (4fn  configuration),  the  stmt  can  be 
cut  off  at  n  3  6  by  the  triangle  rule,  which  governs  the  addition  of 
angular  momenta.  Then  by  observing  which  spherical  tensors  (proportion¬ 
al  to  spherical  harmonics)  are  compatible  with  the  St*  symmetry  of  the 
crystal  field,  we  see  that  the  series  reduces  to  terms 
involving  C®,  C*2,  cP,  C*4,  C*2,  Cp ,  and  C*1*.  Since  the  wave  functions 
have  odd  parity,  parity  considerations  of  the  matrix  elements  eliminate 
the  two  odd  crystal-field  terms,  C*2  and  C*2.  Finally,  the  overall 
phase  can  be  adjusted  so  that  B*1*  is  real.  This  adjustment  reduces  the 
Hamiltonian  to  seven  terms  with  seven  parameters  to  be  determined: 

V_  «  B®C°  +  B®cP  +  B**C**  +  B °<£  +  (rc  b£  +  Im  B^c£  +  xfefc)  .  (2) 

For  Yb3  +  with  a  4f 1 3  configuration,  the  wave  functions  are  for  a 
single  hole.  The  evaluation  of  the  radial  parts  of  the  matrix  elements 
for  Vc  involve  terms  in  rn,  and  these  terms  have  been  calculated  by 
Freeman  and  Watson.^  The  angular  terms  are  the  usual  hydrogenic  func¬ 
tions  •  Intermediate  coupled  states  are  required,  in  general,  for  the 
rare  earth  ions  since  L  and  S  are  not  good  quantum  numbers  due  to  the 
spin-orbit  interaction.  However,  for  Vb3+,  there  is  only  one  term  (2F)j 
thus,  this  consideration  becomes  academic.  The  wave  functions  used  are 
the  same  ones  calculated  by  Pappalardo  and  Wood®  (except  for  a  correc¬ 
tion  of  the  phase  in  the  J  *  5/2  states).  The  14  x  14  perturbation 
matrix  formed  by  these  wave  functions  and  the  Hamiltonian  of  equation 
(2)  can  be  reduced  to  two  identical  4x4  matrices  and  two  identical  3  x 

1 

3  matrices,  the  former  representing  the  four  doubly  degenerate  J  3  7/2 

a 

states  and  the  latter  representing  the  three  doubly  degenerate  J  3  5/2 
states • 


Complete  diagonal izat ion  of  the  contracted  7x7  matrix,  fitting  of 
a  set  of  crystal-field  and  spin-orbit  parameters  to  the  experimental 
data,  and  calculation  of  the  eigenvalues,  wave  functions,  and  g  factors 
were  all  accomplished  by  a  computer  program  written  by  N.  Karayianis  of 
the  Harry  Diamond  laboratories  and  revised  by  this  author,  the  program 
accepted  the  seven  energy  levels  and  the  parallel  g  factors  of  the 
lowest  J  »  5/2  and  J  =*  7/2  states  as  input.  From  these,  it  formed  the 
perturbation  matrix  and  diagonalized  it  by  successive  rotations.  The 
resulting  eigenvalues  were  functions  of  the  six  crystal-field  parameters 
and  X,  the  spin-orbit  parameter,  which  then  had  to  be  fitted  by  a  para¬ 
bolic,  least  squares,  iteration  subroutine.  Goodness  of  fit  was  deter¬ 
mined  by  minimizing  a  quality  factor,  Q,  which  is  a  weighted  root  mean 
square  deviation.  This  quality  factor  is 


Ye*P  -  V 

*i  i 


.calc) 


72 


AX 


exp 


(3) 


0Xp  calc 

where  X^  and  are  the  experimental  and  calculated  values  of  the 
energy  levels  and  g  factors  to  be  compared,  and  Ax®*1*  is  the  experimen¬ 
tal  error  for  these  quantities,  which  is  used  as  a  weighting  factor. 
Where  x  represents  a  g  factor,  AX  contains  an  additional  factor  of  100, 
which  causes  the  Q  for  the  g  factors  to  be  of  the  same  order  of  magni¬ 
tude  as  the  Q  for  the  energy  levels. 

In  many  cases,  the  selection  of  spectral  lines  (energy  levels) 
giving  a  good  fit  was  not  unique.  Equally  good  fits  could  be  obtained 
in  some  cases  by  moving  the  selection  of  the  St«  spectral  line  by  as  much 
as  25  cm"1 .  This  uncertainty  caused  a  problem  in  identifying  the  elec¬ 
tronic  transitions,  but  it  also  led  to  a  deeper  insight  into  the  nature 


of  the  problem.  When  the  matrix  is  diagonalized,  the  resulting  expres¬ 
sions  relating  the  energy  levels  to  the  crystal-field  parameters  are  not 
solvable  analytically.  Even  for  no  J-mixing,  an  attempt  to  solve 
exactly  the  set  of  equations  failed  because  they  were  nonlinear, 
coupled,  and  probably  underdetermined.  The  addition  of  J-mixing  in¬ 
creases  the  underdetermined  nature  of  the  problem.  Thus,  it  is  con¬ 
ceivable  that  the  computer  program  will  calculate  small  values  of  Q  for 
more  than  one  set  of  choices  of  the  seven  St*  lines. 

The  experimental  methods  actually  used  to  identify  the  transitions 
are  discussed  in  I.  The  spectral  lines  and  the  g  factors  tabulated  in  I 
were  entered  into  the  program,  and  the  results  are  shown  in  table  1;  the 
crystal  field  and  the  spin-orbit  parameters  for  the  nine  lattices  will 
reproduce  the  energy  levels  and  g  factors.  As  shown  in  I,  there  were 
some  highly  structured  groups  of  lines  in  the  spectra;  by  selecting 
different  peaks  within  each  group  as  the  S4  line,  the  values  of 
the  B®  changed.  This  variation  in  the  B®  as  the  choice  of  the  S4  lines 
is  varied  among  the  peaks  of  each  group  is  represented  by  the  errors 
given  in  table  1.  No  error  is  given  for  imaginary  (Im)  b£  values  since 
the  variation  was  very  large.  The  values  of  this  parameter  are 
presented  for  completeness  only  and  could  probably  be  set  to  zero  or  any 
other  value  within  several  orders  of  magnitude  with  no  change  in  the 
results . 


TABLE  1.'  FITTED  CRYSTAL-FIELD  PARAMETERS  FOR  NINE  9CHEELITE  CRYSTALS 


b"  CdMoO. 

n  v 

CaW>4 

c*moo4 

Srwo4 

SrHo04 

PbW04 

fbN004 

BaW04 

uyf4 

«7^. 

44«:u 

494*$ 

*-68 

«7-a 

399*2i» 

“5/ 

414*2® 

“25 

399*23 

-35 

*°*-U 

28,^so 

-S20*21 

”u-35 

-340*11 

”U-20 

-311*22 

J  -59 

-556***» 

-20 

K  83,-255 

776^t* 

719*™ 

“115 

739*110 

'”-17 

752*50 

'**“55 

757*55 
•  -8/ 

«i!»; 

559*53 

-35 

n 

32*15 

”-•0 

■S,-M 

-39*15 

-/5 

-uonj 

-106*®° 

-  K  460-w 

452^02 

-155 

552*1  «•/ 
”*-152 

337*0 

■”'-150 

324*05 
32  -103 

269*1 2 / 
-158 

8«0*°® 

-100 

k  *125 

518 

1 

384 

3Y0 

506 

216 

499 

953 

x  -2902* 

-2003* 

-2900* 

-2904* 

-2902* 

-2901* 

-2901* 

-2903* 

-2 

-2897* 

,  not*  1  *  ana  -  valuta  indioat*  arror  limit a. 


This  characteristic  of  Im  b£  was  explained  when  a  series  of  plots 

was  made  showing  the  variation  of  the  energy  levels  when  five  of  the  six 

crystal  parameters  are  held  constant  and  the  sixth  is  varied  by  ±500 

cm'1  about  its  best  fit  value.  Two  of  these  plots  are  shown  as  an 

example  superimposed  in  figure  1.  For  CaWO^ ,  the  three  levels  of  the 

J  =  5/2  manifold  are  plotted  as  a  function  of  B®  and  B®.  The  slope  of 

the  lines,  which  is  a  direct  measure  of  the  dependence  of  the  energy 

levels  on  each  crystal-field  parameter,  is  considerably  greater 

when  B?  is  varied  than  when  bQ  is  varied.  This  difference  illustrates 
Z  b 

the  general  result  in  which  it  was  found  that  the  energy  levels  were 
from  two  to  seven  times  more  dependent  upon  the  B$j ,  B® ,  and  B*|  values 
than  they  were  on  the  sixth  order  terms.  In  other  words,  a  large  fluc¬ 
tuation  in  the  value  of  a  sixth  order  parameter  would  have  a  negligible 
effect  on  the  spectra,  but  a  small  change  in  a  second  or  fourth  order 


ENERGY  (cm1) 


Figure  1.  Dependence  of  CaW04 

J  =  5/2  manifold  on  B®  and  B® . 

Z  b 


VARIATION  IN  CRYSTAL-FIELD  PARAMETER  (cm-1) 
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term  would  cause  large  shifting  in  the  spectral  lines*  Hie  spectrum  is 
almost  completely  independent  of  the  value  of  Zm  .  For  this  reason, 
in  the  following  attempt  to  reproduce  the  crystal-field  parameters  from 
a  theoretical  model,  only  the  second  and  fourth  order  terms  are  con¬ 
sidered,  whereas  the  sixth  order  parameters  are  allowed  to  fall  where 
they  may. 


3.  POINT  CHARGE  MODEL  CALCULATIONS 


Let  us  consider  the  point  charge  model  in  which  the  basic  premise  is 
that  the  impurity  ion  bearing  a  charge  (which  for  the  4f 13  configuration 
of  Yb3+  is  +e)  is  subjected  to  the  electric  field  due  to  the  surrounding 
ions  of  the  lattice.  At  a  point  R  from  the  ion,  this  lattice  charge 
distribution  is  p(R)«  If  the  position  of  the  +e  charge  of  the  ion  is  at 
the  ionic  radius,  r,  then  the  crystal -field  potential  at  that  point  can 
be  expressed  as 


V(r) 


(4) 


/. 


l 

n,m 


4ne 
2n  +  1 


r.  , 

P (R)  —77  V®  (e,,«T)Ym(e,<fr)  dt’ 

n+i  n 


where  (R,8',$’)  and  (r,8,$)  are  the  position  coordinates  of  a  point  in 
the  lattice  charge  distribution,  p(R),  and  the  4fi3  electrons,  respec¬ 
tively;  dr'  indicates  that  the  integration  is  over  th^  lattice  charge 
distribution;  and  r<  is  the  lesser  and  r>  is  the  greater  of  r  and  R.  We 
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will  mium  for  the  present  that  all  of  p (R)  is  external  to  the  impurity 
ion  and  thus  R  >  r .  Also,  we  define  spherical  tensors  to  be  of  the  fora 


Now  we  can  write  the  crystal-field  potential  as 


V(r) 


»,m  LJ  R 


dR'J  rVJ(8,*) 


(5) 


“  ’  l 


AmtrV“ 

n  n 


(6) 


=  I 

n,m 


v“fc“ 

n  n 


(7) 


where  the  V®  are  the  calculated  equivalents  of  the  experimental  crystal- 
n 

field  parameters,  the  .  in  practice,  there  is  an  additional  factor  in 
equation  (6),  so  that  is  written 


(8) 


In  this  equation,  an  is  a  shielding  factor  that  represents  the  extent 
that  the  5s^5p®  electrons  modify  the  field  at  the  4f  electrons  due  to 
the  lattice  ions*  Aa  is  the  electrostatic  potential  due  to  the  ions  of 
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the  crystal  at  the  4f  electrons.  Thus,  the  term  r^d  -  <*n)  pertains  to 

the  impurity  ion,  while  Am  relates  to  the  lattice. 

n 

The  <*n  have  been  calculated  by  Sternheimer  et  al®  for  Pr3+  and  im3+, 
from  which  we  can  obtain  approximate  values  for  Yb3+  by  extrapolation. 
The  second  factor  of  the  is  r11,  which  upon  forming  the  matrix  ele¬ 
ments  becomes  the  expectation  value  ^r11^.  The  calculation  of  the  (rn^ 
for  free  trivalent  rare  earth  ions  has  been  done  by  Freeman  and  Watson7 
using  Hartree-Fock  methods. 

The  lattice  sums  were  calculated  by  changing  the  integral  equation 
(5)  to  a  sum. 


qi  /  (n.  -  m) l 
Rn+1^(n  +  m)  I 


(9) 


This  expression  treats  the  charge  distribution  p(R)  in  equation  (5)  as  a 
set  of  discrete  point  charges,  q^e,  and  sums  over  all  i  of  them  at  their 
positions  (1^,6^,$^).  The  are  spherical  harmonics  as  defined  by  Rose 
or  Edmonds.10'11  Equation  (9)  can  be  reexpressed  in  rectilinear  coordi¬ 
nates  (x,y,z)  with  the  substitution  site  taken  as  the  origin.  The 
coordinates  of  any  ion  (x^,y^,z^)  can  be  expressed  as  functions  of 
the  nearest  neighbor  oxygen  coordinates  (xQ,yQI zQ) .  Using  equation  (9) 
in  equation  (8),  along  with  the  values  of  <rn>  and  (1  -  on) ,  we  calcu¬ 
lated  values  of  with  the  aid  of  a  computer  program  written  by  C.  A. 
Morrision  and  N.  Karayianis  of  the  Harry  Diamond  Laboratories  and  re¬ 
vised  by  this  author.  The  program  performs  the  lattice  sums  of  equation 
(9)  over  all  points  in  the  lattice  for  as  many  "shells,"  that  is,  layers 
of  unit  cells,  surrounding  a  particular  impurity  ion,  as  desired  to 
obtain  convergence. 
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Before  discussing  the  results  of  the  sunning  program,  however,  it  is 

necessary  to  discuss  the  nature  of  this  input  to  the  program.  Beginning 

with  oxygen  coordinates,  it  has  been  pointed  out  in  numerous  places  in 

the  literature  (for  example.  Burns12)  that  small  uncertainties  in  the 

oxygen  coordinates  can  result  in  large  variations  in  the  calculated 

values  of  the  Am .  This  result  was  indeed  borne  out  by  the  calculations 
n 

in  this  work.  Table  2  collects  the  oxygen  (or  fluorine)  coordinates  for 
all  the  crystals  studied  and  gathered  from  various  authors.  For 
example,  in  CaW04  there  are  seemingly  small  discrepancies  between  the 
values  given  by  Kay  et  al13  and  Zalkin  and  Templeton14  and  those  given 
by  Wyckoff.16  However,  they  result  in  a  factor  of  2  difference 
in  and  A £  and  a  factor  of  10  difference  in  A® .  Since  there  is  a  much 
higher  degree  of  accuracy  in  the  neutron  diffraction  studies  of  Kay  et 
al  (which  also  agree  with  the  x-ray  studies  of  Zalkin  and  Templeton), 
the  data  of  Kay  et  al  are  chosen  as  the  preferred  set.  Likewise,  the 
neutron  diffraction  studies  of  Gurmen,  Daniels,  and  King16  on  CaMo04, 
SrW04,  SrMo04,  BaW04 j  of  King  on  LiVF4»  and  of  Leciejewicz1 1  on  PbMo04 
have  been  used  ral'her  than  Hyckoff's  data.  There  are  no  available 
neutron  diffraction  values  of  the  oxygen  coordinates  for  CdMo04  or 
PbWC>4;  thus,  the  sums  for  these;  lattices  were  not  carried  out. 

The  lattice  sums  were  carried  out  by  locating  each  type  of  lattice 
site  with  a  position  vector  as  follows: 

*  ■  *  "J2  *  K’o'Vo)  *  \]! 

(10) 

*  fi[Z(X0'yO'“o)  +  nz]2)l/2  ' 


TABLE  2.  SCHEELITE  COORDINATES 


Crystal 

Rtftrtnet 

*yp*  of 
naasuraasnt 

*o 

*0 

*o 

cawo4 

s 

X-ray 

0.25  t  0.02 

0.11  *  0.02 

0.07  *  0.015 

b 

X-ray 

0.2415  t  0.0014 

0.1504  t  0.0013 

0.006  i  0.0006 

c* 

Mautron 

0.2413  i  0.0005 

0.1511  k  0.0000 

0.006  *  0.0001 

CaMoOj 

a 

X-ray 

0.25  t  0.02 

0.15  t  0.02 

0.07S  t  0.015 

d* 

Mautron 

0.2430  *  0.0010 

0.1459  ±  0.0009 

0.030  t  0.0004 

8rW04 

a 

X-ray 

0.2$  t  0.02 

0.075  t  0.015 

d* 

Mautron 

0.2370  *  0.0000 

0.1107  ±  0.0007 

0.0015  0  0.0003 

8rMo04 

a 

X-ray 

0.25  t  0.02 

0.14  t  0.02 

0.075  *  0.015 

d * 

Nautron 

0.2370  ±  0.0010 

0.1353  t  0.0000 

0.0800  *  0.0004 

PbU04 

a 

X-ray 

0.25  t  0.02 

0.13  ±  0.02 

0.075  t  0.015 

Pb#to04 

a 

X-ray 

0.247  *  0.02 

0.092  *  0.02 

0.0.5  o  o.ots 

a* 

Mautron 

0.2352  ±  0.00060 

0.1134  t  0.00073 

0.0419  *  0.00024 

bswo4 

a 

X-ray 

0.25  t  0.02 

0.11  *  0.02 

0.075  *  O.OtS 

d* 

Mautron 

0.2333  t  0.0005 

0.1214  t  0.0006 

0.0770  *  0.0002 

liyf4 

f 

Mautron 

0.2820  *  0.0011 

0.1642  *  0.0011 

0.0015  4  0.0004 

'lUuu  tuad  in  mu  calmlttiau. 

*»•  ».  G.  Myckoff,  Crystal  Structural,  3,  John  Wlty  and  Son*,  Inc,,  Mem  York  (1965),  19  tt,  taken 
auigly  from  l.  G.  Sillm  and  A,  Hy lander,  Ark.  Kami,  kin.  Gaol.,  IT  A,  do.  4  (19 43). 

°A.  Zalkln  and  0.  H.  Templeton,  J.  (Stem.  Phya.,  #  (1964),  SOI .  aa  quoted  in  o. 

CM.  I.  Kay,  B.  C.  Frazer,  and  I.  Alaadovar ,  J,  Cham,  phy.,  40  (1964),  504. 
dE.  cur  men,  S.  Daniels,  and  J,  S.  King,  3.  Orem.  Phy*.,  55  (1971),  1093-1097. 

®J.  Lsciejeuicz,  Z.  Kriat.,  131  (196S).  1S4-164, 

rJ.  s.  King,  University  of  Michigan,  private  oommuntca t ion . 


where  a  and  c  are  the  lattice  parameters  (given  in  I)  and  X,  Y,  and  Z 

are  the  site  coordinates ,  which  are  expressed  as  functions  of  the 

nearest  neighbor  oxygen  or  fluorine  coordinates  (x0>y0,zo).  The  n  are 

integers  that  run  from  zero  to  some  value  that  is  selected  to  cause  the 

sum  to  converge.  Hie  program  is  written  in  such  a  way  that  the  n  are 

set  at  their  maximum  values  and  stepped  down  in  integral  units  to 

zero.  In  this  way,  the  small  contributions  from  the  outer  lattice  sites 

are  counted  first  so  that  they  do  not  lose  significance  when  added  to 

larger  inner  terms.  In  general,  we  set  the  maximum  vhlue  of  n.  to  be 

half  of  n  and  n  since  c  >  2a.  This  setting  allowed  us  to  carry  out 
x  y 

the  lattice  sums  over  a  roughly  spherical  volume.  Convergence  was 
determined  by  carrying  out  a  number  of  test  sums  for  CaWOj  and  varying 
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the  maximum  values -of  the  (n^n^n^)  to  determine  the  point  beyond  which 
the  lattice  sum  did  not  change.  The  restilts  show  that  the  largest 
contribution  comes  from  the  ligands;  the  ions  farther  away  contribute 
relatively  little.  For  A® ,  a  stable  value  representing  convergence  is 

m 

reached  by  (nx,ny/nz)  -  <o,6,3) ,  while  convergence  for  the  A4 
and  Ag  sums  is  achieved  by  (nx*ny»n2)  *  (4,4,2).  However,  as  a  precau¬ 
tion  against  small  fluctuations— especially  in  A® ,  which  converges  in  an 
oscillatory  fashion,  it  was  decided  to  carry  out  all  sums  to  (n^n^ng) 
=*  (10,10,5). 

In  the  tungstate  and  molybdate  scheelites,  the  bond  between  the 
cation  and  the  heavy  metal  oxide  tetrahedron  is  principally  ionic, 
whereas  the  tetrehedron  itself  is  mainly  covalent.  thus,  the  total 
charge  of  the  tetrahedron  is  -2.0e.  A  molecular  orbital  calculation  by 
Karavelas1 3  on  the  vanadate  tetrahedron  in  a  Cawo^-like  structure  showed 
that  practically  all  the  charge  resided  on  the  four  oxygen  ions;  this 
redistribution  of  charge  leaves  the  vanadium  ion  almost  neutral.  Hie 
extrapolated  ratio  of  charge  for  (WO^)2”  results  in  -0.53e  on  the  oxygen 
ion  and  0.12e  on  the  tungsten  ion.  From  a  suggestion  by  C.  A.  Morrison 
of  the  Harry  Diamond  laboratories,  the  oxygen  charges,  qg,  were  varied 
between  the  limits  -0.5e  and  -2.0e  (the  tungsten  charges  were  varied 
between  the  limits  0.0  and  +6.0)  in  an  attempt  to  find  a  good  calculated 
fit  to  the  experimental  data.  For  U.YF4,  the  (LiF^)3”  tetrahedron  is 
not  covalent,  but  ionic.  Thus,  the  charge  on  the  fluorine  ion  is  not 
expected  to  vary  from  its  value  of  -1.  (For  brevity,  we  refer  to  the 
effective  oxygen  or  fluorine  charge,  qQ,  as  a  dimensionless  number  and 
understand  it  to  be  multiplying  the  electronic  charge,  e.)  In  the 
calculation,  the  values  of  ^r1^  were  from  Freeman  and  Watson,7  and  the 
values  of  an  were  from  Sternheimer  et  al9  and  are  listed  in  table  3. 
The  effective  oxygen  charge  for  CaWO^tYb3*  that  gave  the  best  fit  was 
■  -1.30.  The  crystal-field  parameters  calculated  with  this  charge  and 
also  with  q 0  »  -2.00  are  compared  with  the  experimental  values  in 
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table  4.  The  average  discrepancy  between  the  five  calculated  from 

the  respective  experimental  B  ( ignoring  Im  V?  and  Im  )  was  66  percent 

n  do 

for  ■  -1.30  and  83  percent  for  (Jq  =  -2.00.  The  corresponding  average 
percentage  of  deviation  for  three  parameters  (B® ,  b®,  b^)  is  given  also 
in  table  5. 


TABLE  3.  <r">  AMD  a„  OSEO  IM  CRYSTAL-FIELD 

PARAMETER  CALCULATIONS 


«  <r">‘ 

•»* 

( atomic  units) 

2  0.613 

0.533 

4  0.960 

0.008 

6  3.104 

-0.043 

*A.  J.  Freeman  and  R 

.  E.  Watson, 

Phys.  Rev.,  127  (1962), 

2058. 

bR.  M.  Stemheismr,  . 

M.  Blume, 

and  R.  P.  Pelerla,  Phys 

.  Rev.,  173 

(1968),  376. 


TAIL*  4.  LATTICE  SUMS  OP  DIFFERENT  STUDIES 


< 

CsW04- 

IbMoO^ 

Eraain  at  alc 

This  work 

Sengupta  and  Artaand 

This  «ork 

< 

2880 

4029 

2420 

3916 

< 

-342 

-359 

-200 

-196 

0.86 

-0.52 

-0.3 

-2 

MS  A|| 

-159 

-150 

-47 

-81 

-200 

-237 

-73 

-198 

*•*{ 

-17 

-19 

-4 

-10 

JLL L- 

-15 

-16 

-4 

-10 

Oxygen  coordinates  from  H.  I.  Kay,  8.  C.  Fraser ,  and  I.  Almadovar,  J.  Cham. 
Phys . ,  40_  (1964),  504;  lattice  parameters  from  A.  H.  Morozov  at  al.  Opt. 
Spectrosc.  (USSR),  22  (1967),  139/  oxygen  charge  -  -2.00. 

b Oxygen  coordinates  and  lattice  parameters  from  J.  Lecie jetties,  2.  Krist., 
121  (1965),  158-164 /  oxygen  charge  »  -2.00. 

CM.  V.  Eremin  et  al,  Sov.  Phys.  Solid  State,  11_  (1970),  1697. 

*D.  Sengupta  and  J.  0.  A rtman,  Phys.  B,  1_  (1970),  2986-2988. 
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446 

(-53, 

463) 

-538 

(-59, 

420) 

776 

(-186, 

,  +46) 

-11 

(-38, 

427) 

452 

(-164, 

,  4402) 

518 

Average  deviation  of  fit  of 
B®,  B°,  and  B*  («) 


57 

397 

445  t  10 

-577 

-598 

-587  *  2  , 

815 

586 

581  1  1 

-68 

-6 

i 

-9.1  *  0.6  i 

494 

335 

291  6  1 

65 

-6 

22  t  0.2 

33 

17.8 

11.5 

where  the  sun  is  over  the  three  second  and  fourth  order  terns.  Although 
0  could  be  used  as  a  figure  of  merit,  the  reported  66-  and  83-percent 
fits  were  the  absolute  percentages  of  difference  between  and  b”, 

b“-  V* 

%  difference  *  •■■■■■■• . — 

Bm 
n 

averaged  over  five  parameters,  B°  ,  B*j ,  b£  ,  B°  ,  and  Re  b£  .  It  was  be¬ 
lieved  that,  although  the  routine  to  obtain  a  best  fit  minimized  o,  the 
presentation  of  the  percentage  of  difference  allowed  us  to  visualize  the 
goodness  of  fit  more  easily. 

The  results  of  two  papers  dealing  with  lattice  sums  of  scheelites 
were  examined  for  comparison  with  the  sums  calculated  in  this  work. 
Sengupta  and  Artman13  investigated  neodymium,  Nd3+,  and  neptunium,  TJp'*+, 
in  SrW04,  PbHo04,  and  BaHo04«  Eremin  et  al20  investigated  Nd3+  in 
CaW04,  CdMo04,  SrW04,  PbMo04,  and  BaMo04.  Their  results  are  compared 
with  our  calculations  in  table  4.  We  computed  these  sums  for  comparison 
only,  using  the  lattice  parameters,  the  oxygen  coordinates,  and  the 
oxygen  charges  that  were  used  in  those  two  papers,  and  we  will  show  that 
other  values  for  some  of  these  parameters  are  preferable.  The  monopole 
lattice  sums  calculated  by  Eremin  et  al  agree  well  with  ours,  except 
that  V®  is  only  57  percent  of  our  value.  The  sums  of  Sengupta  and 
Artman  do  not  agree  with  ours,  most  of  ours  being  higher, 
particularly  A^.  Although  neither  paper  details  the  summing  technique, 
the  sums  probably  had  not  converged.  As  is  mentioned  above,  it  was 
found  to  be  necessary  to  sum  from  the  outside  in,  so  that  the  signifi¬ 
cance  in  the  contribution  of  the  outer  term  is  not  l&st  during  trun¬ 
cation  in  the  computer.  It  is  suggested  that  this  method  may  not  have 
been  used  for  those  papers.  Also,  particularly  for  a£,  it  is  necessary 
to  carry  out  the  sum  futher  before  convergence. 


Other  questions  may  be  raised  concerning  the  lattice  subs  calculated 
in  those  two  papers*  Sengupta  and  Artman19  admit  that  the  heavy  metal 
tetrahedra  are  covalently  bonded,  and  thus  its  charge  distribution  is 
altered.  However,  neither  they  nor  Eremin  et  al2f  adjust  the  charges. 
Furthermore,  Eremin  et  al  include  the  dipolar  contributions,  but  only 
from  the  nearest  neighbor  ions.  This  may  not  be  a  valid  treatment  since 
the  dipolar  contribution  from  the  rest  of  the  lattice  may  not  be  negli¬ 
gible.  Also,  Hutchings  and  Ray3  show  that  the  multipolar  series  conver¬ 
ges  slowly.  In  fact,  they  found  the  quadrupole  term  to  be  comparable  to 
the  monopole  term.  Finally,  Eremin  et  al20  used  the  oxygen  coordinate 
data  for  CaW04  to  calculate  the  lattice  sums  for  the  other  four  lattices 
that  they  treated.  In  light  of  the  sensitivity  of  the  sums  to  these 
coordinates,  this  treatment  must  make  their  calculations  for  CdMoo^, 
SrWo4,  FbMoOj,  and  BaMoO^  highly  suspect. 


4.  MODIFICATIONS  TO  CRYSTAL-FIELD  PARAMETER  CALCULATIONS 

It  is  obvious  from  the  foregoing  that  the  point  charge  model,  while 
useful  in  predicting  experimental  results  in  a  qualitative  manner,  is 
not  successful  quantitatively.  Some  suggested  modifications  to  the 
theory  have  been  proposed  including  lattice  polarization3  and  spatial 
extension  of  the  lattice  points2*  and  covalent  bonding23-29  of  the  rare 
earth  ion  in  the  lattice.  These  modifications  have  all  involved  great 
calculational  difficulties  necessitating  approximations  that,  in  turn, 
tend  to  leave  the  modification  suspect. 

Alternatively,  we  have  considered  the  problem  from  the  point  of  view 
that  when  an  Yb3+  ion  enters  a  scheelite  lattice,  certain  physical 
phenomena  can  take  place,  and  an  attempt  has  been  made  to  account  for 
each  effect  in  a  phenomenological  way.  Ms  must  be  wary  of  overparame- 


terizing  the  problem/  but  we  must  still  account  for  as  many  of  the 
physical  phenomena  occurring  as  we  can*  The  introduction  of  these 
phenomenological  corrections  is  not  meant  to  substitute  for  more  exact 
calculations/  but  rather  to  help  determine  which  effects  can  give  agree¬ 
ment  between  calculated  and  experimental  crystal-field  parameters  and  to 
set  limits  on  the  magnitude  of  these  effects. 


*Rie  tungstate  or  molybdate  tetrahedron  is  a  covalent  complex/  and 
the  charge  distribution  on  it  will  be  varied  from  what  it  would  be  if  it 
were  ionic*  Thus,  the  effective  oxygen  charge,  <Jq,  has  been  introduced 
as  an  adjustable  parameter .  (The  charge  on  the  fluorine  ion  ill  LiYF^, 
however,  was  not  varied.)  Varying  <Jq  alone  does  not  provide  a  good 
fit.  Several  other  modifications  have  been  considered  and  are  discussed 
in  the  following  sections. 

4.1  Effect  of  Local  Lattice  Distortion 


A  modification  was  suggested  by  recent  work  on  lattice  distor¬ 
tions  around  impurity  ions  in  alkali  halides  and  alkali  earths.30-32 
Since  the  size  mismatch  between  the  Yb3  +  ions  and  the  divalent  cations 
is  between  8  and  31  percent  in  the  various  scheelite  lattices,  it  is 
expected  that  the  eight  ligand  oxygen  or  fluorine  ions  will  tend  to 
Collapse  somewhat  toward  the  Yb3+  ion.  Thus,  the  value  of  the  lattice 
sums  will  change,  especially  since  the  nearest  neighbor  contribution  has 
been  shown  to  be  the  dominant  term.  There  was  no  a  priori  knowledge  of 
how  far  the  ligands  could  be  allowed  to  move.  It  ms  calculated  that 
there  would  be  sufficient  room  between  the  oxygen  ions  to  allow  the 
ligands  to  collapse  until  they  butted  against  the  Yb3*  ion  similar  to 
hard  sphere  packing.  This  maximum  change  in  ligand  distance  would  range 
between  6  and  16  percent.  Calculating  the  crystal-field  parameters 
using  a  reduced  ligand  distance  “  nRQ,  where  for  CaWO^  n  -  0.942,  we 
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found  a  best  fit  to  occur  for  an  effective  oxygen  charge  of  (Jq  =  1.16. 
The  five-parameter  fit  improved  to  56  percent.  However,  although  a  6— 
to  16-percent  change  in  the  ligand  distance  was  accepted  as 
"reasonable,"  there  was  no  way  of  knowing  that  this  is  the  true  magni¬ 
tude  of  the  distortion.  It  is  probably  safe  to  say  that  some  distortion 
does  occur,  but  it  may  well  be  less  than  the  maximum,  "hard  sphere" 
value  since  the  ligands  sure  part  of  the  covalently  bonded,  heavy  metal 
tetrahedra  and  are  thus  restricted  in  their  motion.  Another  factor  that 
has  not  been  considered  explicitly  is  that,  if  the  ligands  distort, 
there  will  probably  be  distortion,  displacement,  or  rotation  of  the 
near-neighbor  heavy  metal  tetrahedra  to  which  they  are  attached.  Hi  is 
effect  would  change  the  lattice  sums,  but  such  an  effect  would  be 
smaller  than  the  effect  of  moving  just  the  ligands.  Hie  reasons  are 
that  the  contribution  to  the  lattice  sum  from  the  heavy  metal  ion  and 
the  outer  oxygen  ions  on  the  nearest  tetrahedra  is  much  smaller  than 
that  from  the  ligand  oxygen  ions  because  of  their  greater  distance.  The 
amount  of  motion  of  the  outer  ions  will  be  much  less  than  that  of  the 
ligands  because  of  the  binding  of  the  tetrahedra  into  the  rest  of  the 
lattice . 


4.2  Effect  of  Distorton  of  4f  Radial  Wave  Function 


The  next  attempt  at  modifying  the  point  charge  model  proved  to 
be  more  successful  in  providing  a  good  fit  to  the  experimental  data.  We 
considered  the  so-called  "nephelauxetic”  effect,  which  has  been 
suggested  to  occur  when  a  free  ion  is  placed  in  a  crystal.  Because  of 
the  screening  of  the  4f  electrons  by  the  overlapping  charge  clouds  of 
the  surrounding  ligands,  the  4f  electron  wave  function  tends  to  be 
"released"  from  its  nucleus  and  to  expand  slightly  into  the  crystal. 
This  effect  has  been  discussed  in  a  number  of  publications  by 
Jorgensen, 33'3<t  but  there  is  no  quantitative  result  available  that  would 
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indicate  the  amount  of  this  expansion  that  we  might  expect  for  our 
system  of  Yb3+  in  scheelite.  Thus,  we  have  treated  the  radial  wave 
function  as  another  adjustable  parameter  as  follows. 

Freeman  and  Watson^  have  calculated  the  4f  radial  wave  func¬ 
tions  for  the  trivalent  rare  earths  to  be  of  the  form 


-Z./r 

P4£(r)  =  l  C^e  1  ,  i  *  1,  2,  3,  4 


(13) 


with  the  normalization  condition 


P2  (r)  dr  =  1 
4f 


(14) 


For  Yb3+, 


the  values  of  the  and  are  these: 


i 


Ci 


1 

2 

3 

4 


15.287  3914.4363 

8.501  790.99957 

5.667  90.998364 

3.126  4.8064115 


These  yield  the  values  of  given  in  table  3. 
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It  was  apparent  by  observation  that,  if  the  wave  function  did 

indeed  expand,  the  values  of  ^rn^  would  increase  and  the  fit  would 

improve.  The  criterion  for  adjusting  the  wave  function  was  that  it 

should  be  changed  enough  to  shift  the  by  the  desired  amounts,  but 

not  so  much  that  the  probability,  p2  (r),  of  finding  a  4f  electron  at 

4r 

some  distance  from  the  nucleus  would  be  radically  changed.  Having  no 
guidance  on  how  much  the  wave  function  should  expand,  we  sought  merely 
to  keep  it  at  a  minimum.  The  expansion  was  accomplished  by  replacing  z^ 4 
in  equation  (13)  by  Z!^  »  (1  +  K)Z^,  where  K  is  a  constant,  and  by 
some  new  normalization  constant,  «  pC^.  K  and  p  were  chosen  the 

same  for  all  i  =  1,  2,  3,  4.  The  result  of  simultaneously  varying  K  and 
the  oxygen  charge  was  a  notably  improved  fit.  The  best  fit  value  of  K 
and  p  for  CaW04  was  K  =  0.250  and  p  =  0.274.  The  new  wave  function 
parameters  were  these: 


i  C[  -  0.2740^ 


Z'L  =  1.250Zi 


1  4.189 

2  2.329 

3  1.553 

4  0.857 


4893.0454 

988.74946 

113.74796 

6.008014 


The  4f  charge  densities  for  K  =  0,  K=  0.160,  and  K  =  0.250  are  compared 
in  figure  2.  These  values  represent  the  range  over  which  K  varies  for 
all  the  lattices  examined.  The  charge  density  of  a  4f  electron  at  the 
1fb3+  radius  increases  as  K  increases.  The  average  separation  between 
the  calcium  and  oxygen  centers  in  CaW04  is  about  2.45  A.  Thus,  the 
expanded  Yb3+  wave  function  does  not  protrude  too  far  into  the  oxygen 
charge  cloud.  However,  the  new  values  of  ^r1^  for  CaW04  are  greatly 
changed,  those  being 


<r*>  -  1.090,  <r*»>  -  3.035,  and  <r6>  -  17.443 
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RAMUS  (A) 

1.1  1  1.1  t 


Figure  2.  Radial  wave  functions 
of  4f  electrons  of  Yb3+. 

The  best  fit  for  CaWO^  was  obtained  for  sun  effective  oxygen 
charge  of  qg  *  -1.00  with  K  »  0.250.  The  ligand  distance  was  reduced  by 
n  =  0.942  (5.8  percent)  in  this  case. 


The  resulting  average  discrepancy  between  and  (five- 
parameter  fit)  was  17.4  percent.  Considering  only  the  errors 
of  b0  ,  B® ,  and  ,  the  discrepancy  reduces  to  11.5  percent .  Further¬ 
more,  if  the  experimental  error  is  taken  into  account,  the  fit  becomes 
nearly  perfect.  Without  reducing  the  ligand  distance  (n  *  1),  we  can 
get  a  good  fit  also  by  expanding  the  wave  function  somewhat  further  by 
changing  K  from  0.250  to  0.315.  The  three-parameter  fit  (second  and 
fourth  order  terms  only)  is  17.8  percent  without  reducing  the  ligand 
distance . 
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4.3  Some  Other  Possible  Effects 


There  are  yet  two  more  possible  phenomena.  First  it  is  con¬ 
ceivable  that  the  5s25p®  wave  functions  also  would  be  distorted  by  the 
ligands  as  are  the  4f  wave  functions,  this  effect  would  manifest  itself 
by  altering  the  Stemheimer  shielding  parameters,  <*n,  n  =  2,  4,  6.  This 
effect  would  probably  be  small  since  the  5s  and  5p  shells  are  closed, 
whereas  the  4f  shell  is  not.  Furthermore,  any  change  in  the  multipli¬ 
cative  factor  (1  -  otjj)  could  be  taken  up  by  a  change  in  ^r11^,  which  also 

multiplies  Am  to  get  vm. 

n  n 

Second,  no  calculations  have  been  made  to  determine  the  effects 
of  covalency,  overlap,  or  exchange  interactions  between  the  ligands  and 
the  impurity  ions.  If  meaningful  calculations  could  be  made  to  take 
these  effects  into  account,  it  would  not  be  necessary  to  rely  on  best- 
fit  criteria  to  determine  the  extent  of  4f  radial  wave  function 
expression  and  local  lattice  distortion. 


5.  RESULTS  OF  CRYSTAL-FIELD  PARAMETER  CALCULATIONS 

The  lattice  sixn  program  was  run  for  the  seven  crystals  CaW04, 
CaMo04,  SrW04,  SrMo04,  0bMoO4,  BaW04,  and  LiYF4*  The  ion  position  for 
the  other  two  scheelites  studied  in  I  were  not  considered  sufficiently 
accurate  for  detailed  analysis.  Various  values  of  q^,  the  effective 
oxygen  charge,  and  K,  the  radial  wave  function  expansion  constant,  were 
tried.  The  oxygen  and  fluorine  coordinates  used  are  indicated  in  table 
2.  The  program  was  run  also  with  and  without  the  ligand  distance  reduc¬ 
tion  factor,  n*  The  results  are  presented  in  the  following  tables. 
Table  6  gives  the  best  fit  values  of  qQ,  K,  p,  the  resulting  ^rn^,  and 
the  ligand  distance  reduction  factor,  n«  n  was  not  considered  as  an 
adjustable  parameter,  but  only  as  a  multiplicative  factor  of  either 


1*000  for  no  local  distortion  or  sons  value  between  0.835  and  0.943, 
indicating  (for  the  respective  lattice)  complete  relaxation  of  the 
ligands  into  the  impurity  ion's  hard  sphere  radius. 


TABLE  6.  LATTICE  SUM  PITTING  PARAMETERS 


Crystal 

to 

X 

P 

<E»> 

<**> 

<r®> 

q 

Ca«04 

-1.00 

0.250 

0.27401 

1.090 

3.035 

17.443 

0.942 

CaMs04 

-too 

0.240 

0.29084 

1.062 

2.878 

16.110 

0.943 

SrW04 

-1.10 

0.200 

0.36635 

0.958 

2.344 

11.842 

0.893 

Srlto04 

-0.90 

0.230 

0.30846 

0.899 

2.732 

14.895 

0.893 

PbMo04 

-0.90 

0.240 

0.29084 

1.062 

2.878 

16.110 

0.878 

BaW04 

-0.80 

0.250 

0.27401 

1.097 

3.035 

17.442 

0.835 

LiTF4 

-1.00 

0.160 

0.43630 

0.824 

1.929 

8.837 

0.929 

notes: 

qQ  is  the  effective  oxygen  charge.  (For  LiYF^,  it  is  the  fluorine  charge, 
which  does  not  change  from  its  ionic  value.) 

K  is  the  radial  wave  function  expansion  factor. 

p  la  the  normalization  multiplier  for  the  expanded  radial  wave  function. 

q  is  the  ligand  distance  reduction  factor. 

Table  5  gives,  for  CaWO^,  a  comparison  of  the  experimental  crystal- 
field  parameters,  B™,  with  those  calculated  from  the  point  charge 
model  and  its  modifications.  Column  A  gives  the  unmodified  monopole 
lattice  sum  with  q^j  ■  -2.0  and  no  correction  for  ligand  distance  or  wave 
function  expansion  (that  is,  n  *  1.000  and  K  «  0) .  Column  b  is  similar 
to  column  A,  but  the  oxygen  charge  has  been  optimized  to  <Jq  =  -1.30. 
Column  C  has  introduced  the  reduced  ligand  distance  and  has  reoptimized 
the  charge  to  =  -1.16.  Column  D  gives  the  fit  when  qQ  is  held  to 
-0.53  in  accordance  with  the  calculation  of  Karavelas,18  and  the  wave 
function  is  allowed  to  expand.  Column  E  gives  the  fit  ^obtained  with  n  * 
1.000  and  varying  qQ  and  K.  The  best  fit  obtained  is  given  in  column  F, 
which  is  similar  to  the  fit  in  column  E,  but  with  n  “  0.942.  In  columns 
A,  B,  and  C,  the  ^r11^  calculated  by  Freeman  and  Matson7  are  used  and  are 
given  in  table  3.  In  column  F,  the  (rn^  used  are  given  in  table  6.  For 


all  six  columns,  the  shielding  parameters,  aR,  were  extrapolated  from 
Sternheimer  et  al9  and  are  given  in  table  3.  The  progressive  improve¬ 
ment  in  the  fit  is  indicated  by  the  last  row  in  table  5,  where  the 
average  percentage  of  deviation  of  the  fit  of  B^ ,  B9 ,  and  Bjj  is  given* 


Following  this  pattern,  tables  7  to  13  are  similar  to  table  S,  with 
the  results  for  the  other  six  crystals*  Only  two  sums  are  shown  for 
each  crystal  (corresponding  to  columns  A  and  F  in  table  S),  along  with 


the  fitted  B^.  Table  13  gives  the  calculated  odd  lattice  sums  for 
completeness.  They  are  given  for  both  g^  *  -2.0  and  the  optimised  <Jq 
for  each  crystal. 

Considering  three  parameters  (that  is,  eliminating  B9,  Re  B9,  and 

c  m 

Im  B~ ) ,  the  average  difference  between  the  experimental  B^  and  the 

theoretical  for  seven  crystals  is  17  percent.  If  the  experimental 

error  in  the  Bm  is  included,  the  difference  of  fit  reduces  to  8.5 
n 

percent. 


TAB  LI  7.  AGttKBUMT  SMS  CXPIAIMIXTAL  AMD  CALCULATMD  CRYSTAL -FIILO 
PARAMtTIRS  FOR  CBltoO,  LATTtCX 


■ 

n 

0 

494  (-48,  *s> 

0 

-520  (-34,  *21) 

4 

719  (-114,  *79) 

0 

32  (-189,  415) 

111  4 

552  (-152,  +147) 

la  4 

1 

Avmrmgm  deviation  of  fit  of  8® 

8®,  Bj  (%> 


qQ  -  -2.00#  K  -  0  qQ  -  1.10,  K  -  0.240 
1054  533  t  10 

-241  -374  t  2 

240  590  t  1 

-4  -18.2  t  0.6 

68  255  t  1 

_ ^£0 _ 14.3  t  0*2 

54.6  17.3 


lot. I  *  «nd  -  yjIum  indict,  error  limit •  in  B°  eoicn. 


TAB  LI  6*  AGREEMENT  BETWEEN  EXPERIMENTAL  AND  CALCULATED  CRYSTAL -FIELD 
PARAMETERS  FOR  SrWOa  LATTICE 


0 

0 

4 

0 

Re  4 
la  4 


447  ( -68,  +71) 
-340  <-22,  +39) 
739  (-17,  *110) 
-29  (-80,  +14) 
337  (-148,  4 0 ) 
384 


Avtr«9*  deviation  of  fit  of 
B«,  ■«,  and  bJ  («) 


^  *  kn  <*">  (  -  “n) 


2.00,  X  •  0 


1084 

-187 

194 

-7 


qQ  *  -1.10,  K  -  0.200 


457  t  10 
-421  t  2 
489  *  1 
-24  t  0.4 
172  *  1 
35  t  0.2 


17.7 


note:  *  end  -  values  indicate  error  limits  in  tP  column. 


TABLE  9.  AGREEMENT  BETWEEN  EXPERIMENTAL  AND  CALCULATED  CRYSTAL-FIELD 


PARAMETERS  FOR  8rM>04  LATTICE 


n 

g 

B“ 

■>0-1) 

“  -2.00,  It  -  0 

90  •  -0.90.  It  •  0.230 

2 

B 

399  (-57.  *24) 

1058 

374  t  10 

i 

-307  (-54*  +0) 

-167 

-364  t  2 

B 

794  (-91,  *49) 

191 

485  t  1 

I 

-40  (-77,  412> 

-4 

-37  t  0.4 

ft 

273  (-102,  478) 

44 

172  t  1 

9 

I 

310 

-17 

38  t  0.2 

Average  devimtiaa  of  fit  of 

* 

and  bJJ  («> 

70.5 

21.2 

note:  *  end  -  values  indicate  error  lieits  in  tP  col  inn • 

n 
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<■ 

n 

“  2.00.  K  •  0 

^  -  -0.98.  K  -  0.240 

2 

0 

399  (-36.  ♦23) 

1121 

377  ft  10 

4 

0 

-J11  (-.9.  *22) 

-m 

-412  ft  2 

4 

4 

767  (-67.  ♦54) 

187 

539  t  1 

6 

0 

-39  (-76.  ♦  «) 

-6 

-40  ft  0.6 

6 

Ra  4 

324  (-103.  ♦84) 

44 

197  ft  1 

6 

In  4 

216 

-18 

49  ft  0.2 

Aval 

rage  deviation  of  fit  of  1 

>1 

• 

and  b£  (ft) 

75.5 

19.4 

Hotet  ♦  and  -  values  indicate  error  iiaits  in  coluam. 

n 


TAB  LB  11.  AGREEMENT  BETWEEN  EXPERIMENTAL  AND  CALCULATED  CRYSTAL-FIELD 
PARAMETERS  FOR  RaNDf  LATTICE 


n 

n 

* 

2 

0 

404  (-15.  +57) 

4 

0 

-262  (-60,  ♦€) 

4 

4 

661  (-117,  +22) 

6 

0 

-140  (-54,  +17) 

6 

Ra  4 

269  (-168,  T127) 

6 

In  4 

499 

qq  -  -2.00,  K  ■  0  |  qp  »  -6.80,  K  »  0.350 


1057 

-S3 

143 

-7 

26 

-13 


371  ft  10 
-*306  ft  2 
509  ft  1 
-69  ft  0.6 
155  ft  1 
62  ft  0.2 


Average  deviation  of  fit  of 
R®,  8®,  and  Bj  (ft) 


77.1 


16.0 


Vote*  ♦  and  -  values  indicate  error  Units  in  B*  colunn. 

n 
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TAB  LI  12.  AGREEMENT  BETWEEN  EXPERIMENTAL  AMD  CALCULATED  CRYSTAL-FIELD 
PARAMETERS  FOR  LiYFj  LATTICE 


0 

0 

4 

0 

Re  4 
X*  4 


281 

<-50. 

♦40) 

-556 

<-20, 

♦34) 

569 

(-36, 

♦53) 

-106 

<-«7. 

♦200) 

840 

(-100 

♦200) 

953 

Average  deviation  of  fit  of 
B»,  B«,  and  Bj  <»> 


-1.00,  K  -  0 


-133 

232 

-1 

52 

-8 


48*5 


go  »  -1*00,  K 


0.160 


284  t  10 
-432  t  2 
585  t  1 
-11  t  0.6 
283  t  1 
35  t  0.2 


14.6 


Note i  +  and  -  values  indicate  error  limits  in  0*  column. 


TUU  13.  CALC01AT1D  VRLOZS  FOR  ODD  LKRICK  BIM 


cryatal 

<*0 

ft*  k\ 

0.  *| 

Rk  a| 

CaWO. 

-2.00 

224 

310 

-146 

16 

-1.00 

-45 

-234 

-93 

-8 

CJUtoO. 

-2.00 

204 

254 

-143 

14 

-1.10 

32 

-199 

-102 

-8 

SrWO. 

-2.00 

204 

242 

-117 

14 

-1.10 

-10 

-316 

-102 

-12 

SrltoO. 

-2.09 

317 

205 

-115 

13 

-0.90 

-0 

-301 

-08 

-13 

fbttoO. 

-2.00 

320 

236 

-113 

14 

-0.90 

-3 

-323 

-91 

-14 

Bam. 

-2.00 

201 

157 

-07 

11 

•0.00 

-■ 

-349 

-03 

-17 

utf4 

-1.00 

107 

-46 

-171 

1 

Die  greatest  error,  in  general,  lies  in  the  n  »  6  terms.  The  reason 
for  this  error  may  be  partly  that  <^r^  is  considerably  more  sensitive  to 
perturbation  by  the  ligand  ions  than  is  ^r2^  or  ^r*1^  •  It  is  more 
sensitive  because  the  sixth  power  amplifies  the  contribution  of  ^r11^  of 
the  radial  wave  function  tail,  which  extends  further  into  the  lattice 
than  do  the  tails  of  the  second  and  fourth  powers. 

6.  CONCLUSIONS 

The  point  charge  model  in  its  basic  form  does  not  accurately  repro¬ 
duce  the  experimental  data.  However,  modifying  the  theory  by  including 
various  effects  that  occur  when  an  impurity  ion  enters  a  lattice  yields 
greatly  improved  results.  In  particular,  accounting  for  the  covalent 
bonding  of  the  (W04)2“  and  (Mo04)2“  tetrahedra  by  adjusting  the 
effective  oxygen  charge  and  accounting  for  the  expansion  of  the  4f 
radial  wave  function  ( nephelauxetic  effect)  seem  to  be  especially  useful 

in  fitting  the  theoretical  V*®  to  the  experimental  B® .  The  results 

n  n 

strongly  suggest  that  4f  wave  functions  expand  radially  when  the  ion  is 
introduced  into  a  solid.  Local  distortion  of  the  ligand  ions  may  or  may 
not  occur.  A  good  fit  between  calculated  and  experimental  crystal-field 
parameters  is  obtained  whether  or  not  we  assume  local  distortion. 
However,  neglect  of  this  effect  requires  the  assumption  of  a  greater  4f 
wave  function  distortion  to  obtain  an  equally  good  fit  to  the 
experimental  crystal -field  parameters. 

Furthermore,  adjusting  the  oxygen  charge,  considering  the  local 
lattice  distortion  at  the  ligands,  and  examining  the  nephelauxetic 
effect  are  all  manifestations  of  the  covalent  nature  of  the  lattice  and 
the  bonding  of  the  impurity  ion  in  the  lattice.  Nevertheless,  they  have 
all  been  treated  in  a  phenomenological  manner  and  are  thus  only  an 
approximate  approach  to  the  problem.  However,  rigorous  covalency  calcu- 


lations  are  usually  quite  difficult  and  require  assumptions  to  be  made 
as  to  the  positions  and  the  wave  functions  of  the  ions  concerned.  They 
are  often  no  more  accurate  than  the  phenomenological  treatment  of  this 
work.  The  results  of  this  study  also  confirmed  the  fact  that  the  posi¬ 
tions  of  the  ions  in  the  lattice  bear  a  critical  relation  to  the  lattice 
sums  from  which  the  crystal-field  parameters  are  calculated. 
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